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We study thermoelectric effects in superconducting nanobridges and demonstrate that the mag- 
nitude of these effects can be comparable or even larger than that for a macroscopic circuit. It is 
shown that a large gradient of the electron temperature can be realistically created on nanoscale 
and masking effects of spurious magnetic fields are minimal in nanostructures. For these reasons 
nanodevices are favorable for studying the thermoelectric effect in superconductors. 
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I. INTRODUCTION 

The discrepancy between the theory and experiment 
concerning the thermoelectric phenomena is a long stand- 
ing problem in physics of superconductors. The ther- 
moelectric phenomena in the superconducting state were 
first discussed by Ginzburgi as early as 1944. In the 
presence of a temperature gradient, there appears in a 
superconductor a normal current of the form given by 



j T = -i)VT 



(1.1) 



where rj is the corresponding transport coefficient. In the 
bulk of a homogeneous isotropic superconductor the total 
electric current must vanish, and as was pointed out by 
Ginzburgi, the normal current is offset by a supercurrent 
j s so that the total current in the bulk 



jT+j S = 0. 



(1.2) 



This makes impossible standard studies of the thermo- 
electric effect in a homogeneous isotropic superconductor. 
Ginzburg considered also simply-connected anisotropic 
or inhomogeneous superconductors as systems where it 
is possible to observe thermoelectric phenomena by mea- 
suring the magnetic field generated by a temperature gra- 
dient. 

Theory of the effect was further developed in 1970s 2 . It 
was noted in particular that the offset supercurrent is re- 
lated to a difference of the order parameter phases within 
the simple-connected superconductor. This phase differ- 
ence can be measured either in superconducting interfer- 
ometer or in the loop formed by different superconduc- 
tors where a magnetic flux is generated in the presence 
of a temperature difference. It stimulated experimental 
study of the effect. Although the first experiment per- 
formed by Zavaritsky^ is in a rather good agreement with 
the existing theory, further experiments (see e.g4£) ex- 
hibit temperature-dependent magnetic fluxes five order 
of magnitude larger than is predicted by the theory^. A 
possibility to observe large thermoelectric fluxes is dis- 
cussed in£ and is related to the phonon drag effect near 
the contact of the two superconductors with different val- 
ues of superconducting gap. However, the predicted en- 



hancement factor, the ratio of the Fermi energy and De- 
bye energy, is not big enough to bridge the gap between 
the experiments, and the theory. 

From the experimental point, the main difficulty is due 
to the fact that the thermoelectric effect is small at low 
temperatures, and one needs to single it out from vari- 
ous masking effects. The most obvious one is related to 
the temperature dependence of the magnetic field pene- 
tration length 7-8 . As a result, in the presence of a back- 
ground magnetic field, the magnetic field within the su- 
perconductor is temperature dependent. This can mask 
the genuine thermoelectric effect. It is important to note 
that later on it was shown that the co-existence of a tem- 
perature gradient and a supercurrent leads to variation of 
the gauge invariant scalar potential (f> related to an imbal- 
ance between the electron-like and hole-like quasiparticlc 
branches2ii2iiiii2. In contrast to the thermoelectric flux, 
the experimental studies of this effect were in agreement 
with the theory^. 

The goal of the present paper is to discuss the geome- 
try of experiment where the thermoelectric effect is par- 
ticularly large while the masking effects are significantly 
suppressed. Therefore one can hope that it will guar- 
antee unambiguous measurements of the thermoelectric 
effect. With this purpose in mind we will consider ther- 
moelectric effect in superconducting circuit containing a 
point contact. 

Thermoelectric phenomena in superconducting nanos- 
tructures have some unique specific features that (i) are 
favorable from the experimental point of view and (ii) re- 
quire certain revision of the existing theory. Experimen- 
tally, the advantage is that one is able to create very large 
temperature gradients so that the intrinsic thermoelectric 
current becomes larger and easier to observe. We note 
such a favorable possibility can be realized only in the 
systems where the electrons can be heated as compared 
to the lattice. Indeed, a realistic nanostructures imply a 
presence of an (insulating) substrate it is deposited on, 
and thus large gradients of the lattice temperature can 
not be achieved because of the phonon heat conductance 
in the substrate. In contrast, the electron temperature, 
according to the Wiedemann-Franz law, behaves in the 
same way as the electrostatic potential. A detailed dis- 
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cussion demonstrating possibility to obtain large gradi- 
ents of electron temperature in metal nanostructures is 
given in the Appendices. Most important, the parasitic 
effects due to the trapped magnetic field are much less 
pronounced in small size structures. 

On the theoretical side, a revision is needed because 
the earlier theories have considered bulk samples. Their 
sizes have been assumed to be much larger than the char- 
acteristic lengths such as the London penetration length 
and the length at which the offset supercurrent is gener- 
ated. When applied to bulk samples, there is no need to 
specify and go into detail of the mechanism by which the 
normal thermoelectric current is converted into the offset 
supercurrent. This approach is valid provided the sample 
size much exceeds the size of the region where the normal 
thermoelectric current is converted into the offset super- 
current, ft is well known from the microscopic theory 
that the conversion occurs as a result of branch-mixing 
scattering processes, where electron-like excitations are 
scattered to the hole branch of the excitation spectrum 
and vice versa. Microscopic mechanism of the branch- 
mixing is known to be inelastic scattering, impurity scat- 
tering in the case of anisotropic superconductors, and 
Andreev reflection if inhomogeneity of the order parame- 
ter gap is present. If the bulk scattering is the dominant 
mechanism, the conversion takes place along the branch- 
mixing diffusion length Lb, Lb = \ADt&, D and r& being 
the diffusion constant and the branch-mixing relaxation 
time, respectively. For a nanostructure of the size com- 
parable with the branch-mixing length, the standard the- 
ory of the thermoelectric phenomena (that assumes local 
compensation of the thermoelectric current) is not ap- 
plicable. Indeed, in this case the normal thermoelectric 
current can be offset also by a normal diffusion current 
rather than by a supercurren1jA£. 

In addition to kinetics, there are important differences 
in electrodynamics of superconducting nanostructures. 
In particular, it is related to the so-called kinetic induc- 
tance £k 



£ k = ^#- (1-3) 



Here Al is the London penetration length, L is the cir- 
cuit length while S is the circuit cross-section. is 
inversely proportional to S and may be larger than the 
magnetic inductance of the thermoelectric loop for very 
small values of S. In this case, the local compensation 
of the current in Eq. (|I.2|) turns out to be energetically 
unfavorable and the electrodynamical part of the theory 
requires a revision too. 

In what follows we will develop a theory of thermo- 
electric effect in superconducting nanostructures. It will 
include the above kinetic and electrodynamical effects. 



II. CHARGE IMBALANCE DISTRIBUTION 

Consider two superconducting films (banks) connected 
by a narrow wire of the length L and cross-section S; 
the transverse sizes of the wire are assumed to be much 
smaller than the London penetration length A^. In this 
case, the current is distributed homogeneously across the 
wire crossection, and the problem is one dimensional. De- 
note x the coordinate along the wire and choose the ori- 
gin in the middle of the wire. We analyse a diffusive wire 
and assume that the temperature varies linearly along the 
wire between, its values at the banks being Tl and Tr. 
Note that the thermoelectric current is considered a con- 
stant equal to — 17VT in the wire and zero in the banks. 
This assumption holds for 3D structures where both the 
temperature gradient and electric current density quickly 
decays within the contact. Naturally, we assume that the 
thickness of the wire is much smaller than the thicknesses 
of the banks. 

First, we briefly overview the well-known physics of 
branch imbalance in superconductors. The total electric 
current density, j = j s +j n , is a sum of the supercurrent, 
j s , and normal, j n , components. The supercurrent reads 



where p s is the superfluity momentum, 

P S = ?V*--A, (2.2) 
2 c 

X and A being the phase of the order parameter and the 
vector potential, respectively. 
The normal current, 

jn=3T+3D, (2.3) 

is a sum of the thermoelectric current jV = — f?VT, and 
the diffusion component, jo = — erV</>, related to the 
branch imbalance specified by the gauge invariant poten- 
tial <j) as , 

0=^X + ^, (2.4) 
2e 

if being the scalar potential. In the vicinity of the critical 
temperature, the diffusion current is proportional to the 
normal state conductance a. 

A. Diffusion limit 

The potentials p s and (j) are found from the continu- 
ity equation divj = 0, and the equation which describes 
transformation of the normal current into supercurrent 
that results in the following equation for (f> in the wire 
(see, for instance, Ref. 

V 2 0-^ = O (2.5) 
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where 77, is the branch imbalance relaxation time while 
Lb = y/Drb is the branch imbalance relaxation length^. 
If the banks are made of superconductors with different 
gap values, Eq. 1)2. 5f) requires a boundary condition^ to 
account for the Andreev reflection at the interface. The 
latter plays the role of a surface mechanism of branch- 
mixing. 

The boundary condition to Eq. (|2.5|) rather generally 
takes the formic 



1 

a 



x=±L/2 



1 

± 5 



(2.6) 



x=±L/2 



where j n is the a;— component of the normal charge cur- 
rent in Eq. Q2.3[l . and Cb is an effective relaxation length 
controlled by the Andreev reflection at the wire-bank in- 
terface as well as the branch-mixing rate in the banks. 

Solution to Eq. (|2.5|l with the boundary condition 
Eq. (213) reads 



4>(x) = (f) b 



sinh 7rh 
sinh 757^ 



In the present paper, we are interested in the opposite 
limit of a short wire, L « Lb- It follows from Eq. i|2.7|) 
that the potential difference in this limit is 



Jt 



1 



A<p = 

For a short enough wire, L < Cb, we obtain 
Acj) ~ —AT 



(2.10) 



(2.11) 



In this case of a short superconducting wire, the thermo- 
electric potential difference A<fi is of the order of that in 
the normal state. 

In a short wire, the supercurrent Eq. I|2.8|l is homoge- 
nous, 



L 



3 s 



JO- JT 



L 



(2.12) 



where as before jo is the total electric current through 
the wire. 



where 



JT 



rr + srrcoth^ 



(2.7) 



The condensate momentum p s and the supercurrent 
j s is found from the continuity equation divj = 0, that 
is j(x) = jo, where jo is a constant. From the condi- 
tion j s (x) +j n (x) = jo where j n (x) = j T - aVcp(x), the 
distribution of the supercurrent is given by the following 
expression 



js{x) - jo = - jT 1 - 



1 



cosh 



2L b 



^ + coth ^- 

JLb *L, b 



sinh 



2L b 



(2. 



If there is no electrical connection between the banks 
other than the wire, the total current must be zero, 
jo = 0. Otherwise, the constant jo is found from electro- 
dynamical considerations considered in the next section 
Thus the temperature difference between the banks 
leads to a creation of a potential difference equal to 
Acf) — 2<j>b with 4>b in Eq. (|2.7() . The potential difference 
can be measured, for instance, as described in Ref. UtI 
Note that this effect if of the nature considered by Arte- 
menko and VolkoAii^. However, they treated a macro- 
scopic circuit with a size L much larger than Cb so that 
the potential difference was concentrated near the inter- 
face region thus involving only a small part of the total 
temperature difference AT = LVT. As a result, their 
estimate for "thermoelectric" potential difference is 



Acj) ~ — AT— 
a L 



(2.9) 



for a macroscopic wire, the length of which L exceeds the 
microscopic scales Lb and Cb- 



B. Ballistic bridge 

When we studied the branch imbalance in the previous 
section, for the simplicity we have exploited the diffusive 
approximation. However, the largest values of jr cor- 
respond to the largest values of the mean quasiparticle 
free path within the wire (leading to larger rf). So one 
expects the largest effect for a ballistic bridge. In this 
case one can estimate the quasiparticle thermoelectric 
current with the help of a procedure similar to the one 
used in Ref. 18. Namely, one has in mind that the quasi- 
particle distribution function within the ballistic bridge 
is formed by the quasiparticles entering the bridge from 
the banks. One also notes that the distribution function 
is constant along the quasiparticle trajectory while the 
quasi-equilibrium distribution functions of the left and 
right banks correspond to different temperatures (T^ and 
Tr, respectively). Thus for the quasiparticle distribution 
function, Fb, within the bridge one has 



Fn 



e v 



f{t l ) + e 



■v x - I F{T R ). (2.13) 



Here v x is the component of the electron velocity along 
the bridge direction, and F(Tl,r) stands for the equilib- 
rium distribution function corresponding to the temper- 
ature Tl^r. We have taken into account that the (group) 
quasiparticle velocity differs from the "bare" electron ve- 



locity by a factor £/e = £ P /y Cp + A 2 , £ p being the ki- 
netic energy counted from the Fermi energy, and A being 
the superconductor energy gap. Given the distribution 
function in Eq. I|2.13|l . the normal thermoelectric current 
reads 



JT 



= e^v x F E 
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As usual in the theory of thermoelectric phenomena, the 
contributions of electrons and holes to the current tend to 
cancel each other, and the net effect is sensitive to details 
of the band structure and to the energy dependence of the 
density of states, in particular. At temperatures T > A, 
the order of magnitude of the thermoelectric current can 
be estimated as 



j T ~ ev F n 



T 2 R ) 



(2.14) 



where vp and ep are the Fermi velocity and energy, re- 
spectively, and n is the electron density. For a rough 
estimate, assume that the temperature difference is com- 
parable to T c . In this case, 



jt ~ env F 



(2.15) 



Thus a presence of a temperature drop at the contact 
between two superconducting banks leads to formation 
of the thermoelectric current through the nanobridge the 
order of magnitude of which can be evaluated according 
to Eq. H2.15(l . The total thermoelectric current, It, is 
found by the multiplication of the current density jx and 
the bridge cross-section S, It — JtS. 



III. THERMOELECTRIC FLUX 

As we have discussed above, we study a nanostructure 
that consists of a superconducting bridge with a thick- 
ness and a width smaller than the London penetration 
depth X L . The bridge joins two banks made of the same 
superconductor (with a critical temperature T cl and a 
thickness smaller than Al). By means of a point- like 
heating, using e.g. N-S tunnel junction (see Appendix), 
the banks are kept at different temperatures. We assume 
that the bridge region carrying thermoelectric current It 
is short-circuited by superconducting branch with sizes 
larger than Al forming a loop of the linear size C. The 
behaviour of the system is different for the two limiting 
cases: a) C > L^, b) C < Lb. 

We start our analysis with the first one, that is the 
case when the branch imbalance relaxation length Lb is 
much shorter than the size of the system. It can be re- 
alized in particular if the near-contact region at least for 
one of the banks is covered by the superconductor with a 
larger gap leading to effective imbalance relaxation due 
to Andreev reflections. If the circuit is simple-connected 
the thermoelectric current is compensated by the super- 
current created due to the Andreev reflection or bulk 
mechanisms of the charge imbalance relaxation. Thus 
the decay of the normal thermoelectric current is locally 
compensated by supercurrent. 

The situation becomes different if the circuit is not 
simple-connected, i.e., when another branch (made, for 
instance, of the material with a larger T c ) closes the 
loop. In this case, the net electric current, built of the 



normal thermoelectric and superconducting components, 
through the bridge may be finite for the charge current 
continuity is maintained by the supercurrent I c through 
the branch closing the loop: I c is actually the electric 
current circulating in the loop, and (I c — It) is the su- 
perconducting component of the net current through the 
bridge. The circulating current I c can be readily evalu- 
ated minimizing the total energy W of the system. The 
latter is given by the following expression, 



W= \{It-Ic 



) 2 C k + \llL 



(3.1) 



The first term originates from the kinetic energy of su- 
perconducting electrons in the bridge, £/- being the well- 
known kinetic inductance, 



C k 



LA| 



(where L and S, as above, are the bridge length and 
cross-section, respectively). The second term in Eq. (|3.1|) 
is the energy of magnetic field created by the circulating 
current I c , and C is the inductance of the loop, which is 
close to the geometrical inductance of the macroscopic 
branch. Minimizing W with respect to I c , one obtains 



c k + c 



and, thus the thermoelectric magnetic flux is 

C k C 



$r = h 



c k + c 



(3.2) 



(3.3) 



The flux <I>t is controlled by the smaller of the induc- 
tances in question. 

Note that if C k <C C, $t does not depend on C and 
is estimated as 



I T C k 



LXj 
S 



(3.4) 



In the dirty limit, the penetration depth Al is related to 
its value in thebulk pure material as 

where £o ~ Vp/A is the coherence length, and l e is the 
electron elastic mean free path. As it can be seen, this 
result coincides with the predictions of the papers^ as- 
suming that the thermoelectric current is almost com- 
pletely compensated by the supercurrent. 

However, the situation is qualitatively different if 
Ck 3> C, a condition which can be realistically met 
for a nanoscale bridge. Indeed, assuming L ~ y/S and 
L ~ l e ~ 10 _6 cm, for Ao ~ 10~ 5 cm, £o ~ 10~ 4 cm, one 
obtains Ck ~ 10 -2 cm. This means that the kinetic in- 
ductance Ck may be comparable to the magnetic geo- 
metric inductance C even for a relatively large, nearly 
macroscopic loop. In this case, the normal thermoelec- 
tric current generated by the bridge is non-locally short- 
circuited by the supercurrent through the macroscopic 
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branch than being offset locally by the supercurrent. In 
this limit, one has from Eq. 1)3. 3[) : 



d> 7 



I T C. 



(3.5) 



Despite the absence of the current cancellation within the 
bridge, the flux through the loop is Ck/ C times smaller 
than it has been predicted in earlier papers^. At the 
same time, the magnetic field within the structure prac- 
tically coincides with its value for a normal metal struc- 
ture. Correspondingly, it can be much larger than for the 
thermoelectric effect in macroscopic circuits. Indeed, as- 
suming that the inductance C is of the order of the linear 
size of the circuit, our estimates for the magnetic field 
from Eqs. (j3.4jl . and l|3.5|l are 



Hi 



Ck 
'£2 



and 



c k ^>c 



(3.6) 



(3.7) 



Thus the " thermoelectric" magnetic field is the larger the 
smaller is C and is much larger for the regime Ck > C 
than for a "macroscopic" considered earlier in24. We 
believe that this factor significantly suppresses a possible 
role of masking effects. 

In the limiting case (b), when the charge imbalance 
length Lf, is much shorter than the size of the system 
(C < Lb), the quasiparticle thermoelectric current is 
not converted into a supercurrent but short-circuited by 
the normal current through the closing branch (as it oc- 
curs in normal metal thermoelectric circuits). The nor- 
mal charge current in the loop generates a magnetic flux 
which in turn generates a circulating supercurrent I c in 
the direction opposite to the normal current. In this case, 
the energy W = I 2 C k /2 + (I T - I c ) 2 C/2 is built of the 
supercurrent kinetic energy I 2 Ck/2 and the magnetic en- 
ergy (It — I c ) 2 C/2. Minimizing W, 



Ir = h 



c 



C + C k ' 



and the total thermoelectric flux, <E>t = {It ~ Ic)C, is 
again given by Eq. (|3.3|l . Therefore, the thermoelectric 
flux $t is completely controlled by the normal compo- 
nent if Ck > C. 

Let us estimate the largest possible values of $t which 
can be realized for large C. We have 



l 2l 

S 



(3.8) 



where L and 5* are the bridge length and cross-section 
respectively. 

Correspondingly, 



®t ~ envF — 



/' 



le 



(3.9) 



In what follows we will assume that all the sizes of 
the bridge are of the same order while one should also 
put l e ~ L. Assuming T/(i ~ 1CT 4 (T ~ IK), A 2 ~ 
l(r 10 cm 2 , £ ~ 10~ 4 cm one has $ T - l(r 3 $ . 

For smaller C the magnetic fluxes are smaller than the 
above estimate but the magnetic fields are higher. 



IV. CONCLUSION 

In this paper we have analyzed the thermoelectric ef- 
fects in superconducting nanostructures. When the size 
of a thermoelectric circuit is less than the branch imbal- 
ance length, the very picture of the thermoelectric effects 
becomes different from that considered earlier for macro- 
scopic systems: rather than being offset locally by the 
supercurrent, the quasiparticle thermoelectric current is 
short-circuited nonlocally, by the diffusion current in the 
branch closing the circuit, similar to the picture of the 
effect in normal metal thermoelectric loops. The thermo- 
electric effects in superconducting nanostructures may be 
comparable with that in systems of macroscopic size sys- 
tems. At the same time, the masking effects inherent for 
macroscopic superconductors can be eliminated so that 
nanoscale structures are promising for studying the ther- 
moelectric effects in superconductors. 
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APPENDIX A: TEMPERATURE DISTRIBUTION 

The purpose of this section is to discuss the conditions 
when one can ascribe different temperatures to elections 
in two banks connected by a short bridge. 

In recent years it has been demonstrated that the elec- 
tronic temperature of a metal film may substantially dif- 
fer from the lattice temperature of the dielectric sub- 
strate. For quasi-2D metallic nanostructures at low tem- 
peratures, there are two factors that are favorable for 
such a possibility. First, small electron-phonon collision 
rates prevent effective transfer of heat to the phonon sys- 
tem of the substrate. Second, the phonon heat conduc- 
tivity of the substrate at small spatial scales turns out 
to be smaller than the electron heat conductivity within 
the films since the phonon mean free path is limited by 
the spatial inhomogeneity. Using the Wiedemann - Franz 
law, one estimates the electron heat current within the 
metal layer of a length L and cross-section S as 



Qe 



AT 



-SD e h 



(Al) 
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where AT is the temperature difference, D e is electron 
diffusivity, and n ~ (T/ep)n is the concentration of 
quasiparticles participating in the heat transfer (where n 
is the total electron concentration while ep is the Fermi 
energy). At the same time, the heat flux from the film 
to the substrate can be estimated as 



^sub 



SLhAT 

7~e— ph 



(A2) 



where T c _ p h is electron-phonon relaxation time. From 
Eqs. HAljl. and (|A2|) . one sees that Q e \ > Q su b provided 

1? < l e V F T c -ph- 

It is also instructive to compare the electronic heat flux 
Q e i with the heat flux Q ph supported by phonons in the 
substrate and "shunting" the electron flux. One easily 
obtains that Q e i > Q p h provided 



L s mm(w, l p u) 

d VFle 



T 
Tn~ 



El 

T 



< 1 



(A3) 



where To is the Debye temperature of the substrate, w 
is the width of the metal layer, d is the layer thickness, s 
is the sound velocity while l e and l p h are the mean free 
paths of electrons within the layer and phonons within 
the substrate, respectively. Since the electron heat con- 
ductivity dominates provided any of the aforementioned 
conditions holds, one concludes that at low temperatures 
the electron temperature is mainly controlled by electron 
heat conductivity of the metal structure. 

Consider now a point ballistic bridge connecting to 
metal banks with different electronic temperatures. It 
follows from the above considerations that the tempera- 
ture drop is concentrated mainly within the contact re- 
gion. Indeed, for 3D geometry and a diffusive transport 
in the bulk, the temperature distribution in the banks 
near the bridge follows the same law as an electric po- 
tential distribution, that is the temperature drop is con- 
centrated in the bridge. If the whole structure is made of 
a metal film of the same thickness and with the diffusive 
electron transport this statement holds only with a loga- 
rithmic accuracy because of the 2D character of electron 
diffusion. However if the thickness of the bridge region is 
much smaller than the thicknesses of the banks (that is if 
the configuration is a 3D-like one) the temperature drop 
is again completely restricted by the contact region. The 
same holds provided the electron transport within the 
contact and near-contact regions is ballistic. Indeed, it 
follows from the fact that under the Wiedemann-Franz 
law the temperature profile is similar to the electric po- 
tential profile while in 2D ballistic structures the poten- 
tial drop is concentrated in the contact region. 

It is expected that very large values of AT can be 
realized in the point contact geometry. Indeed, one can 
apply for the heat flux the same arguments as for electric 
current through the point contact, namely, that the 
relaxation processes for the electrons take place within 
the bulk of the sample at distances (~ DeT^) 1 / 2 . Thus 
enormous values of temperature gradient and heat flux 
density do not lead to destruction of the bridge. 



1. Electron heating 

Let us consider the important practical question con- 
cerning the generation of the temperature gradient across 
the bridge. We have assumed above that the excitations 
within the one of the banks are heated as compared to 
the excitations in another one. Since we deal with super- 
conductors, it excludes the Joule heating. On the other 
hand, microwave heating implies relatively large areas. 
In our opinion, the best way is to heat electrons on one 
on the banks using a tunnel S-I-N junction. The junc- 
tion is formed by a normal film of area S2 put on the 
top of the superconducting bank (with a thin insulating 
layer). When the bias eV across the S-I-N junction is 
much larger than the superconductor energy gap, high- 
energy electrons tunnelling from N layer will relax mainly 
due to creation of electron-hole pairs within the supercon- 
ducting layer. To have the electron temperature formed, 
one should have 

S 2 > V F l e Tee, 

T ee being the electron-electron scattering time. 

Now let us compare the thermal current from the 
heated superconducting layer to the substrate and the 
thermal current through the point contact to the "cold" 
bank. Assuming that the thickness of the superconduct- 
ing layer forming the tunnel junction and the layer form- 
ing the point contact are the same, one finds that the 
thermal current through the contact dominates provided 



S 2 L/w < leVFT^-pfr 



(A4) 



where L and w are the point contact length and width, 
respectively. If L ~ w one concludes, that this condition 
can hold since at low enough temperatures T ep h > r ee . 
Correspondingly, in this case only a region with the area 
S2 (under the tunnel junction) is heated with respect to 
the rest of the device, the heat leak being due to ther- 
mal current through the point contact. Certainly, one 
should also assume that the area of the superconducting 
layer in the "cold" bank is large enough to ensure effi- 
cient heat withdrawal to its substrate. In this case one 
easily obtains 



AT ~ IV- 



wdD„nT 



(A5) 



where / is the current through the tunnel junction. 

The main conclusion following from the considerations 
given above is that it is possible to have "point-like" elec- 
tron heating restricted by the area ~ VFleT ee - Its linear 
dimensions for realistic estimates can be as small as 3 
/im. Correspondingly, if the inductance loop has macro- 
scopic size this local heating (and corresponding local 
variation of the penetration length) is not expected to af- 
fect the temperature-dependent (or rather l^-dependent) 
flux through the loop. 
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